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Abstract 

The dynamics of typical phase transitions is studied out of equilibrium in weakly 
coupled inflaton-type scalar field theories in Minkowski space. The shortcomings of 
the effective potential and equilibrium descriptions are pointed out. A case of a rapid 
supercooling from Tj > T c to Tf <C T c is considered. The equations of motion up 
to one-loop for the order parameter are obtained and integrated for the case of "slow 
rollover initial conditions" . It is shown that the instabilities responsible for the process 
of phase separation introduce dramatic corrections to the evolution. Domain formation 
and growth (spinodal decomposition) is studied in a non-perturbative self-consistent 
approximation. For very weakly coupled theories domains grow for a long time, their 
final size is several times the zero temperature correlation length. For strongly coupled 
theories the final size of the domains is comparable to the zero temperature correlation 
length and the transition proceeds faster. We also obtain the evolution equations for 
the order parameter and the fluctuations to one-loop order and in a non-perturbative 
Hartree approximation in spatially flat FRW cosmologies. The renormalization, and 
leading behavior of the high temperature limit are analyzed. 



1 Introduction and Motivation 



Inflationary cosmological models provide very appealing scenarios to describe the early evo- 
lution of our universe [fj]]. Since the original model proposed by Guthp|, several alternative 
scenarios have been proposed to overcome some of the difficulties with the original proposal. 

Among them, the new inflationary model[|3], |], ^, ^] is perhaps one of the most attractive. 
The essential ingredient in the new inflationary model is a scalar field (the "inflaton") that 
undergoes a second order phase transition from a high temperature symmetric phase to a 
low temperature broken symmetry phase. The expectation value (or thermal average) of the 
scalar field serves as the order parameter. Initially at high temperatures, the scalar field 

1 Based on Lectures Delivered at the XIII Course of the International School of Cosmology and Gravitation, 
Erice, May 1993, and at the 17th Johns Hopkins Workshop on Current Problems in Particle Theory, July 
1993. 
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is assumed to be in thermal equilibrium and 0^0. The usual field- theoretic tool to study 
the phase transition is the effective potential [0, j|, |§. 

At high temperatures, the global minimum of the effective potential is at (f) = 0, whereas 
at low temperatures there are two degenerate minima. 

The conjectured behavior of the phase transition in the new inflationary model is the 
following: as the universe cools down, the expectation value of the scalar field remains close 
to zero until the temperature becomes smaller than the critical temperature at which the 
effective potential develops degenerate minima away from the origin. When this happens, 
the scalar field begins to "roll down the potential hill". In the new inflationary scenario, 
the effective potential below the critical temperature is extremely flat near the maximum, 
and the scalar field remains near the origin, i.e. the false vacuum, for a very long time and 
eventually rolls down the hill very slowly, hence the name "slow rollover" . During this stage, 
the energy density of the universe is dominated by the constant energy density of the false 
vacuum V e ff(4> = 0), and the universe evolves rapidly into a de Sitter space (see for example 
the reviews by Kolb and Turner fTOfl, Linde|ll[j and Brandenberger [|12|) . Perhaps the most 
remarkable consequence of the new inflationary scenario and the slow rollover transition is 
that they provide a calculational framework for the prediction of density fluctuations [jl3[| . 
The coupling constant in the typical zero temperature potentials must be fine tuned to a 



very small value to reproduce the observed limits on density fluctuations [0, |TT . 

The slow rollover scenario is based on the static effective potential. Its usefulness in a time 
dependent situation has been questioned by Mazenko, Unruh and WaldjjTJJ who argued that 
the dynamics of the cooling down process is very similar to the process of phase separation 
in statistical mechanics, and that the system will form domains within which the scalar field 
will relax to the values at the minima of the potential very quickly. It is now accepted that 
this picture may be correct for very strongly coupled theories but not for weakly coupled 
scenarios as required in inflation. 

Guth and Pi[l5|], studied the effects of quantum fluctuations on the time evolution below 
the critical temperature by treating the potential near the origin as an inverted harmonic os- 
cillator. They recognized that the instabilities associated with these upside-down oscillators 
lead to an exponential growth of the quantum fluctuations at long times and to a classical 
description of the probability distribution function. Guth and Pi also recognized that the 
static effective potential is not appropriate to describe the dynamics, that must be treated 
time dependent process. 

Subsequently, Weinberg and Wu [|l6[1 , have studied the effective potential, particularly in 
the situation when the tree level potential allows for broken symmetry ground states. In this 
case, the effective potential becomes complex. These authors identified the imaginary part 
of the effective potential with the decay rate of a particular initial state. 

In statistical mechanics, the imaginary part is a consequense of a sequence of thermody- 
namically unstable states. 

The effective potential offers no information on the dynamics of the process of the phase 
transition. As mentioned above it becomes complex in a region in field space corresponding 
to thermodynamically unstable states. 
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Furthermore, in an expanding gravitational background the notion of a static effective 
potential is at best an approximate one. To understand whether one may treat the problem 
in the approximation of local thermodynamic equilibrium, the time scales involved must be 
understood carefully. In a typical FRW cosmology the important time scale is determined 
by the Hubble expansion factor H(t) = a(t)/a(t) whereas the equilibration processes are 
determined by the typical collisional relaxation rates T(T) ~ X 2 T in typical scalar theories. 
Local thermodynamic equilibrium will prevail if T(T) ^> H(t). In de Sitter evolution, after 
a typical phase transition at T « T c « 10 14 Gev, H « 10 _5 T c (ll| and local thermodynamic 
equilibrium will not prevail for weakly coupled inflaton theories in which A ~ 10~ 12 — 10~ 14 [|TT1 
11. 

Thus in these scenarios, weakly coupled scalar field theories cannot be assumed to be in 
local thermal equilibrium. Typical expansion rates are much larger than typical equilibration 
rates and the phase transition will occur very rapidly. The long-wavelength fluctuations 
will be strongly out of equilibrium as they typically have very slow dynamics (see[T7|] and 
references therein). The phase transition must be studied away from equilibrium. 

A complete discussion of these issues and the complex effective potential may be found 
in the articles by Boyanovsky and de Vega|18| and Boyanovsky et. al.fT/ 



2 Non-Equilibrium time evolution 

Let us consider the situation in which for time t < the system is in equilibrium at an 
initial temperature Tj > T c where T c is the critical temperature. At time t — the system 
is rapidly "quenched" to a final temperature below the critical temperature Tf < T c and 
evolves thereafter out of equilibrium. 

What we have in mind in this situation, is a cosmological scenario with a period of rapid 
inflation in which the temperature drops very fast compared to typical relaxation times of 
the scalar field. 

Precisely because of the weak couplings and critical slowing down of long-wavelength 
fluctuations, we conjecture that an inflationary period at temperatures near the critical 
temperature, may be described in this "quenched" approximation. Another situation that 
may be described by this approximation is that of a scalar field again at Tj > T c suddenly 
coupled to a "heat bath" at a much lower temperature (below the transition temperature) 
and evolving out of equilibrium. 

The "heat bath" could consist of other fields at a different temperature. 

To understand whether this quenched approximation is valid or not will require a deeper 
understanding of the initial conditions. 

Although we are currently studying the case of inflationary cosmologies, we will here 
concentrate on the dynamics of a supercooled phase transition in Minkowski space. 

This situation is modelled by introducing a Hamiltonian with a time dependent mass 
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term 



Hit) = |^ 3 x|in 2 (x) + i(V<l>(a;)) 2 + im 2 (t)<|. 2 (a;) + A$4 (x) | (1) 

m 2 (t) = m 2 9(-t) + (-/i 2 )0(t) (2) 

where both m 2 and \j? are positive. We assume that for all times t < there is thermal 
equilibrium at temperature T», and the system is described by the density matrix 

Pi = (3) 
Hi = H(t<0) (4) 

In the Schroedinger picture, the density matrix evolves in time as 

p{t) = U(t)p l U-\t) (5) 

with U(t) the time evolution operator. 

An alternative and equally valid interpretation (and the one that we like best) is that the 
initial condition being considered here is that of a system in equilibrium in the symmetric 
phase which is then evolved in time with a Hamiltonian that allows for broken symmetry 
ground states, i.e. the Hamiltonian ([l], |^) for t > 0. 

The expectation value of any operator is thus 

< O > (t) = Tre-^U- x (t)OU(t) /Tre~ PiHi (6) 

This expression may be written in a more illuminating form by choosing an arbitrary time 
T < for which U(T) = exp[— iTHi] then we may write exp[— /^i/j] = exp[— iHi(T — ifa — 
T)] = U(T — z/3j,T). Inserting in the trace f/ _1 (T)[/(T) = 1, commuting £7 -1 (T) with 
and using the composition property of the evolution operator, we may write (§) as 

<0>(t)= TrU(T - ifii, t)OU(t, T)/TrU(T - ifc, T) (7) 

The numerator of the above expression has a simple meaning: start at time T < 0, evolve to 
time t, insert the operator O and evolve backwards in time from t to T < 0, and along the 
negative imaginary axis from T to T — ifa. The denominator, just evolves along the negative 
imaginary axis from T to T — ifa. The contour in the numerator may be extended to an 
arbitrary large positive time T' by inserting U(t,T')U(T' ,t) = 1 to the left of O in (0), thus 
becoming 

<0>(t) = TrU{T - ifr, T)U{T, T')U{T', t)OU(t, T)/TrU(T - if3 u T) (8) 

The numerator now represents the process of evolving from T < to t, inserting the operator 
O, evolving further to T", and backwards from T' to T and down the negative imaginary 
axis to T — ifli. Eventually we take T — > — oo ; T' — > oo. It is straightforward to generalize 
to real time correlation functions of Heisenberg picture operators. 
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As usual, the insertion of an operator may be achieved by inserting sources in the time 
evolution operators, defining the generating functionals and eventually taking functional 
derivatives with respect to these sources. Notice that we have three evolution operators, 
from T to T', from T', back to T (inverse operator) and from T to T — i(3i. 

Since each of these operators has interactions and we want to generate the diagrammatics 
of perturbation theory from the generating functionals, we use three different sources. A 
source J + for the evolution T — > T', J~ for the branch T' — > T and finally J@ for T — > T—ify. 
The denominator may be obtained from the numerator by setting J + = J~ = 0. Finally the 
generating functional Z[J + , J", J 13 } = TrU(T - i/3 h T; jf>)U(T, V; J~)U(T\ T; J+), may be 
written in term of path integrals as (here we neglect the spatial arguments to avoid cluttering 
of notation) 



Z[J + , J', .f] = J DQD^Dfy J V^ + V^'V¥e^r {^ + ^-^-^\} 



X 



(9) 



with the boundary conditions $+(T) = $P(T-ipi) = $ ; $+(T') = $ _ (T') = $ 2 ; $~(T) = 
$^(T) = $i. As usual the path integrals over the quadratic forms may be done and one 
obtains the final result for the partition function 

Z[J + J' J 13 } e^ ^ ' dt [ Cm ^- iS/5J+) ~ Cmt ^ &/5J ' ) ]\ e {* fr~ tfh dtc ™t(- i5 / SJI3 )} x 

e {Uc dtl J c rf *2J c (tl)Jc(t 2 )G c (tl,t 2 )} 

Here J c are the currents defined on the segments of the contour [ 18] J ± , J^[^| and G c is 



the Green's function on the contour (see below), and again the spatial arguments have been 
suppressed. 



In the two contour integrals (on t\] t 2 ) in ( P~0| ) there are altogether nine terms, corre- 
sponding to the combination of currents in each of the three branches. However, in the limit 
T — > — oo, the contributions arising from the terms in which one current is on the (+) or 
(— ) branch and another on the imaginary time segment (from T to T — ifii), go to zero when 
computing correlation functions in which the external legs are at finite real time. For theses 
real time correlation functions there is no contribution from the J 13 terms. These cancel 
between numerator and denominator, and the information on finite temperature is encoded 
in the boundary conditions on the Green's functions (see below). Then for the calculation 
irrelation 

z[j + ,r 



of finite real time correlation functions the generating functional simplifies to[p4|, [25[ 

Jyi dt[c in t(-i6/&J+)-C in t(iS/8J-)]} ^ 
Jy^Sr dtl St ' dt 2 Ja(ti)J b (t2)G ab (t 1 ,t 2 )'j 



with a, b = +, — . 

This formulation in terms of time evolution along a contour in complex time has been 
used many times in non-equilibrium statistical mechanics. To our knowledge the first to 
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use this formulation were Schwinger|26|] and Keldysh|27j (see also Mills[f28j). There are 
many articles in the literature using these techniques to study time dependent problems, 



some of the more clear articles are by Jordan [31], Niemi and Semenoff |23| , Landsman and 
van Weert [ 29] , Semenoff and Weiss |30[|, Kobes and Kowalski||32||, Calzetta and Hu |[25| ,Paz[|33 



and references therein (for more details see [15, |TT|, |T9 



The Green's functions that enter in the integrals along the contours in (|H], [TT]) are given 
by (see above references) 



G ++ (t 1} t 2 ) 

G~(ti,t 2 ) 
G + -(h,t 2 ) 

G<(T,t 2 ) = G>(T-iPi,t 2 ) 



G>{h, t 2 )e(h - 1 2 ) + G < (t 1 , t 2 )e{t 2 - h) 



G > (t 1 , t 2 )Q(t 2 - h) + G < (t h t 2 )e(t! - t 2 ) 
-G<(t u t 2 ) 

-G > (t 1 ,t 2 ) = -G < (t 2 ,t 1 ) 



(12) 
(13) 
(14) 
(15) 
(16) 



As usual G < ,G > are homogeneous solutions of the quadratic form with appropriate 
boundary conditions. We will construct them explicitly later. The condition fllED is recog- 



nized as the periodicity condition in imaginary time (KMS condition) ||34||. 
To obtain the evolution equations we use the tadpole method0, and write 

$±(£,t) = 000 + ^ ± (f,t) 



(17) 



Where, again, the ± refer to the branches for forward and backward time propagation. The 
reason for shifting both (±) fields by the same classical configuration, is that <p enters in the 
time evolution operator as a background c-number variable, and time evolution forward and 
backwards are now considered in this background. 

The evolution equations are obtained with the tadpole method by expanding the La- 
grangian around <fi(t) and considering the linear, cubic, quartic, and higher order terms in 
\l/ ± as perturbations and requiring that 

< ^ ± (x,t) >= 0. 

It is a straightforward exercise to see that this is equivalent to extremizing the one-loop 
effective action in which the determinant (in the logdet) incorporates the boundary condition 
of equilibrium at time t < at the initial temperature. 

To one loop we find the equation of motion 



d 3 k 
W) 3 



(-i)G k (t,t) = 



with Gk(t,t) = G^(t,t) = G^(t,t) is the spatial Fourier transform of the equal-time Green's 
function. 

Notice that 



'-iG k (t,t)) =< 9t{t)*+M)> 
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is a positive definite quantity (because the field ^ is real) and as we argued before (and 
will be seen explicitly shortly) this Green's function grows in time because of the instabilities 
associated with the phase transition and domain growth [| 



IE] 



These Green's functions are constructed out of the homogeneous solutions to the operator 
of quadratic fluctuations 



£ + S» + 



Ul 



M 2 (t) 



A 



(m 2 + -tfM-t) 







The boundary conditions on the homogeneous solutions are 



u£(t < o) 

w<(fc) 



A 



k 2 + m 2 + -</>■ 
2 



(19) 
(20) 

(21) 
(22) 



where fa is the value of the classical field at time t < and is the initial boundary condition 
on the equation of motion. Truly speaking, starting in a fully symmetric phase will force 
fa = 0, and the time evolution will maintain this value. Therefore we admit a small explicit 
symmetry breaking field in the initial density matrix to allow for a small fa. The introduc- 
tion of this initial condition seems artificial since we are studying the situation of cooling 
down from the symmetric phase. However, we recognize that the phase transition from the 
symmetric phase occurs via formation of domains (in the case of a discrete symmetry) in- 
side which the order parameter acquires non-zero values. The domains will have the same 
probability for either value of the field and the volume average of the field will remain zero. 
These domains will grow in time, this is the phenomenon of phase separation and spinodal 
decomposition familiar in condensed matter physics. Our evolution equations presumably 
will apply to the coarse grained average of the scalar field inside each of these domains. 
This average will only depend on time. Thus, we interpret (fi as corresponding to the coarse 
grained average of the field in each of these domains. The question of initial conditions on 
the scalar field is also present (but usually overlooked) in the slow-rollover scenarios but as 
we will see later, it plays a fundamental role in the description of the evolution. 

The identification of the initial value tpi with the average of the field in each domain is 
certainly a plausibility argument to justify an initially small asymmetry in the scalar field 
which is necesary for the further evolution of the field, and is consistent with the usual 
assumption within the slow rollover scenario. 

The boundary conditions on the mode functions U k (i) correspond to "vacuum" boundary 
conditions of positive and negative frequency modes (particles and antiparticles) for t < 0. 

Finite temperature enters through the periodicity conditions (pl|) and the Green's func- 
tions are 



G>(t,t') 
G<(t,t>) 



2w<(£;) 1 - e 
G > (t',t) 



-f3 i u) < (k) 



(23) 
(24) 
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The effective equations of motion to one loop that determine the time evolution of the 
scalar field are 



d 2 (f)(t) 
dt 2 



+ m 2 {t)<p(t) + 



dt 2 



+ k 2 + M 2 (t) 



Ul 







d 3 k U£{t)U k {t) 
(2tt) 3 2w<(A;) 



coth 



(25) 
(26) 



with (p , (|2J. 

This set of equations is too complicated to attempt an analytic solution. They will be 
dealt with numerically. However, before doing this, we should note that there are several 
features of this set of equations that reveal the basic physical aspects of the dynamics of the 
scalar field. 

i) : The effective evolution equations are real. The mode functions U k {t) are complex 
conjugates of each other as may be seen from the time reversal symmetry of the equations, 
and the boundary conditions (pT|). This situation must be contrasted with the expression 
for the effective potential for the analytically continued modes. 

ii) : Consider the situation in which the initial configuration of the classical field is near 
the origin fa rs 0, for t > 0. The modes for which k 2 < (k max ) 2 ; (k max ) 2 = fi 2 — -fa? are 
unstable. 

In particular, for early times (t > 0), when <f. 
imately as 



0, these unstable modes behave approx- 



Then the early time behavior of 

1 



A k e Wkt + B k e 

(Z4 + (t))* 



w k t 



1 - i 



W, 



B, = At 



2 ^ j3 



A- 2 



(27) 
(28) 

(29) 
(30) 



:-iG k (t,t)) 



2cj<(A;) 



1 



iG k (t,t)) is given by 

cosh(2W k t) - 1] 



/i 2 + m 2 



^ ~ f # - k 2 



coth[(3 i u < (k)/2] (31) 



This early time behavior coincides with the Green's function of Guth and Pi[]15] and Weinberg 
and Wu||16|| for the inverted harmonic oscillators when our initial state (density matrix) is 
taken into account. 

Our evolution equations, however, permit us to go beyond the early time behavior and 
to incorporate the non-linearities that will eventually shut off the instabilities. 
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These early-stage instabilities and subsequent growth of fluctuations and correlations, are 
the hallmark of the process of phase separation, and precisely the instabilities that trigger 
the phase transition. 

It is clear from the above equations of evolution, that the description in terms of inverted 
oscillators will only be valid at very early times. At such times, the stable modes for which 
k 2 > (k max ) 2 are obtained from (|27|) , (|^) , ( p9| ) by the analytic continuation Wh — > 

-iu;>(Jfe) = p-/i 2 + |0^. 

For t < 0, U^{t)U^{t) = 1 and one obtains the usual result for the evolution equation 

d^jt) | dv eff {4>) _ Q 

dt 2 d(f) 

with V e ff(4>) the finite temperature effective potential but for t < there are no unstable 
modes. 

It becomes clear, however, that for t > there are no static solutions to the evolution 
equations for (fi(t) ^ 0. 

iii) Coarsening: as the classical expectation value <fr(t) "rolls down" the potential hill, 
4> 2 {t) increases and (k max (t)) 2 = jj 2 — ^4> 2 {t) decreases, and only the very long-wavelength 
modes remain unstable, until for a particular time t s ; {k max (t s )) 2 = 0. This occurs when 
4> 2 {t s ) = 2/z 2 / A, which is the inflexion point of the tree level potential. In statistical mechanics 



this point is known as the "classical spinodal point" and t s as the "spinodal time" [21], [22 
When the classical field reaches the spinodal point, all instabilities shut- off. From this point 
on, the dynamics is oscillatory and this period is identified with the "reheating" stage in 
cosmo logical scenarios [11], |T^| . 

It is clear from the above equations of evolution, that the description in terms of inverted 
oscillators will only be valid at small positive times, as eventually the unstable growth will 
shut-off. 

The value of the spinodal time depends on the initial conditions of 4>(t). If the initial value 
(pi is very near the classical spinodal point, t s will be relatively small and there will not be 
enough time for the unstable modes to grow too much. In this case, the one-loop corrections 
for small coupling constant will remain perturbatively small. On the other hand, however, 
if (pi « 0, and the initial velocity is small, it will take a very long time to reach the classical 
spinodal point. In this case the unstable modes may grow dramatically making the 
one-loop corrections non-negligible even for small coupling. These initial conditions 
of small initial field and velocity are precisely the "slow rollover" conditions that are of 
interest in cosmological scenarios of "new inflation" . 

The renormalization aspects have been studied in reference [18H and we refer the reader 
to that article for details. 
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3 Analysis of the Evolution 



As mentioned previously within the context of coarsening, when the initial value of the scalar 
field <pi ~ 0, and the initial temporal derivative is small, the scalar field slowly rolls down the 
potential hill. But during the time while the scalar field remains smaller than the "spinodal" 
value, the unstable modes grow and the one-loop contribution grows as a consequence. For 
a "slow rollover" condition, the field remains very small (0 2 (t) <C 2/i 2 /A) for a long time, 
and during this time the unstable modes grow exponentially. After renormalization, the 
stable modes give an oscillatory contribution which is bound in time, and for weak coupling 
remains perturbatively small at all times. 

Then for a "slow rollover" situation and for weak coupling, only the unstable modes 
will yield to an important contribution to the one-loop correction. Thus, in the evolution 
equation for the scalar field, we will keep only the integral over the unstable modes in the 
one loop correction. 

Phenomenologically the coupling constant in these models is bound by the spectrum of 
density fluctuations to be within the range X R ~ 10~ 12 — 10 _14 [|TTJ. The stable modes will 



always give a perturbative contribution, whereas the unstable modes grow exponentially in 
time thus raising the possibility of giving a non- negligible contribution. 

With the purpose of numerical analysis of the effective equations of motion, it proves 
convenient to introduce the following dimensionless variables 



t = ^ R t;q = k/fi R 
n \i) = ; L 2 = + fog 

V^R ^R 



(32) 
(33) 



To account for the change from the initial temperature to the final temperature (Tj > 



T c ; Tf < T c ) we parametrize [35 



m 



MO) 

Mo) 



T 2 

- 1 

c 

1 L 

7^2 

c 



(34) 
(35) 



where the subscripts (R) stand for renormalized quantities, and —^r(0) is the renormalized 
zero temperature "negative mass squared" and T 2 = 24/i|,(0) / Xr. Furthermore, because 
(^max(^)) 2 — \i\ and Ti > T c , for the unstable modes % > (k max (t)) and we can take 
the high temperature limit coth[j3 i Lj < (k) /2] w 2T i /u < (k). Finally the effective equations of 
evolution for t > 0, become, after using oj\ = ^ 2 R {(f + L 2 ) and keeping only the unstable 
modes as explained above (g 2 < (5 mra (r)) 2 ), 



dr 2 



t)(t)-v(t) + rf{r) + 
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(wM) M+(t)U-(t) 

2 9 V 7 q y ' dq = 
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g 2 + L 2 



ci 2 

^ + g 2 -(g m ax(t)) 2 



For Tj > T c and 1/ <C T c the coupling 



l"3r/ 2 (r) 



v / 6Ar Ti 



2tt 2 



T 2 

1 - -i- 

1 T? 



(36) 

(37) 
(38) 
(39) 



i is bound within the range g ~ 10 -7 — 10~ 8 . The 
dependence of the coupling with the temperature reflects the fact that at higher temperatures 
the fluctuations are enhanced. 

From fl5E|) we see that the quantum corrections act as a positive dynamical renormalization 
of the "negative mass" term that drives the rolling down dynamics. It is then clear, that the 
quantum corrections tend to slow down the evolution. 

In particular, if the initial value r](0) is very small, the unstable modes grow for a long 
time before t](t) reaches the spinodal point i](r s ) = l/y/3 at which point the instabilities 
shut off. If this is the case, the quantum corrections introduce substantial modifications to 
the classical equations of motion, thus becoming non- perturbative. If 77(0) is closer to the 
classical spinodal point, the unstable modes do not have time to grow dramatically and the 
quantum corrections are perturbatively small. 

Thus we conclude that the initial conditions on the field determine whether or not the 
quantum corrections are perturbatively small. 

Figures (1,2) depict the solutions for the classical (solid lines) and quantum (dashed lines) 
evolution. 

For the numerical integration we have chosen L 2 = 1, the results are only weakly de- 
pendent on L, and taken g = 10~ 7 , we have varied the initial condition 77(0) but used 
Mill — n 

We recall, from a previous discussion that 77 (r) should be identified with the average of 
the field within a domain. We are considering the situation in which this average is very 
small, according to the usual slow-rollover hypothesis, and for which the instabilities are 
stronger. 

In figure (1) we plot rj vs r for g = 1(T 7 , 77(0) = 2.3 x 1(T 5 r/(0) = ; L = 1. The solid 
line is the classical evolution, the dashed line is the evolution from the one-loop corrected 
equation of motion. We begin to see that the quantum corrections become important at 
t ~ 10/ fin and slow down the dynamics. By the time that the classical evolution reaches 
the minimum of the classical potential at 77 = 1, the quantum evolution has just reached 
the classical spinodal point r] = l/y/3. The quantum correction becomes large enough to 
change the sign of the "mass term" |18[| , the field continues its evolution towards larger values, 
however, because the velocity is different from zero. As 77 gets closer to the classical spinodal 
point, the unstability shuts off and the quantum correction arising from the unstable modes 
become small. From the spinodal point onwards, the field evolves towards the minimum and 
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begins to oscillate around it. The quantum correction will be perturbatively small, as all the 
instabilities had shut-off. Higher order corrections, will introduce a damping term as quanta 
may decay into elementary excitations of the true vacuum. 

Figure (2) shows a dramatic behavior for 77(0) = 2.258 x 1CT 5 ; = for the same 
values of the parameters as for figure (1). The unstable modes have enough time to grow 
so dramatically that the quantum correction becomes extremely large overwhelming the 
"negative mass" term near the origin. The dynamical time dependent potential, now becomes 
a minimum at the origin and the quantum evolution begins to oscillate near rj = 0. The 
contribution of the unstable modes has become non-perturbative, and certainly our one- loop 
approximation breaks down. 

As the initial value of the field gets closer to zero, the unstable modes grow for a very long 
time. At this point, we realize, however, that this picture cannot be complete. To see this 
more clearly, consider the case in which the initial state or density matrix corresponds exactly 
to the symmetric case, rj = is necessarily, by symmetry, a fixed point of the equations of 
motion. Beginning from the symmetric state, the field will always remain at the origin and 
though there will be strong quantum and thermal fluctuations, these are symmetric and will 
sample field configurations with opposite values of the field with equal probability. 

In this situation, and according to the picture presented above, one would then expect 
that the unstable modes will grow indefinetely because the scalar field does not roll down 
and will never reach the classical spinodal point thus shutting-off the instabilities. What is 
missing in this picture and the resulting equations of motion is a self-consistent treatment 
of the unstable fluctuations, which must necessarily go beyond one- loop approximation. 



4 Domain formation and growth: 

The instabilities correspond to the growth of long- wavelength fluctuations and the formation 
of domains within which the field is correlated. These domains will grow as long as the 
instabilities persist. In order to understand the process of domain growth when the average 
value of the field remains zero it proves illuminating to understand the tree level correlations. 
The relevant quantity of interest is the equal time correlation function 

S(f;t) = ($(f,t)$(0,t)> (40) 
S(f;t) = J -0j- 3 e iUr S(k-t) (41) 

S(k;t) = (*z(t)*_ u (t)) = (-iGF(t;t)) (42) 

where we have performed the Fourier transform in the spatial coordinates (there still is 
spatial translational and rotational invariance). Notice that at equal times, all the Green's 
functions are equal, and we may compute any of them. 

Clearly in an equilibrium situation this equal time correlation function will be time 
independent, and will only measure the static correlations. In the present case, however, 
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there is a non trivial time evolution arising from the departure from equilibrium of the 
initial state. 

The unstable contribution to the Green function at equal times is given by ([H]), in this 
case for fa = 0. 

It is convenient to introduce the following dimensionless quantities 



K 



k_ 

m f 



m 2 __ 


[Tf 


- T 2 ] 




c 





t = mft ; x = rriff 



(43) 



Furthermore for the unstable modes k 2 < mj, and for initial temperatures larger than the 
critical temperature T 2 = 24/i 2 /A, we can approximate coth[/3ju;<(/c)/2] w 2T i /uj < (k). Then, 
at tree-level, the contribution of the unstable modes to the subtracted structure factor (|42| 
SW(h,t) - S^(k,0) = (l/m f )S^(K,r) becomes 



<S ( °W) 



W K 




k 2 + L 2 



-f, [cosh(2iy K r) - 1] 



(44) 

(45) 
(46) 



To obtain a better idea of the growth of correlations, it is convenient to introduce the 
scaled correlation function 



V(x,t) 



A 



6m 2 



% s k 2 dk sin(fcr) 
o 2tt 2 (kr) 



[S(k,t)-S(k,0)} 



(47) 



The reason for this is that the minimum of the tree level potential occurs at A$ 2 /6m 2 = 1, 
and the inflexion (spinodal) point, at A$ 2 /2m 2 = 1, so that T>(0, r) measures the excursion 
of the fluctuations to the spinodal point and beyond as the correlations grow in time. 

At large r (large times), the product k 2 S(k, t) in ( i7|) has a very sharp peak at K s — 1/ \fr. 
In the region x < y/r the integral may be done by the saddle point approximation and we 
obtain for T//T c ^ the large time behavior 



V(x,r) « P ( 0,r ) e X p[--] 1 l^i 



P(0,r) 



A 



12tt 3 



'ft) 3 | exp[2r] 



T 2 



1] 



3 

r 2 



(48) 
(49) 



Restoring dimensions, and recalling that the zero temperature correlation length is £(0) = 
l/y/2/ji, we find that for Tj »s the amplitude of the fluctuation inside a "domain" ($ 2 (t)), 
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and the "size" of a domain £d(£) grow as 



m)) w exp[^(0)] 

« (8V2)^(0)X^ (51) 

An important time scale corresponds to the time r s at which the fluctuations of the 
field sample beyond the spinodal point. Roughly speaking when this happens, the insta- 



bilities should shut-off as the mean square root fluctuation of the field y($ 2 (t)) is now 
probing the stable region. This will be seen explicitly below when we study the evolution 
non-perturbatively in the Hartree approximation and the fluctuations are incorporated self- 
consistently in the evolution equations. In zeroth order we estimate this time from the 
condition 3V(0,t) = 1, we use A = 1CT 12 ; Tj/T c = 2, as representative parameters (this 
value of the initial temperature does not have any particular physical meaning and was 
chosen only as representative). We find 

t s « 10.15 (52) 
or in units of the zero temperature correlation length 

t w 14.2f (0) (53) 



for other values of the parameter t s is found from the above condition on ([49]). 

Clearly any perturbative expansion will fail because the propagators will contain unstable 



wavelengths, and the one loop term will grow faster than the zero order, etc.|T7|]. We now 
turn to a non-perturbative analysis (for details see [[l7| ). 

As the correlations and fluctuations grow, field configurations start sampling the stable 
region beyond the spinodal point. This will result in a slow down in the growth of cor- 
relations, and eventually the unstable growth will shut-off. When this happens, the state 
may be described by correlated domains with equal probabibility for both phases inside the 
domains. The expectation value of the field in this configuration will be zero, but inside each 
domain, the field will acquire a value very close to the value in equilibrium at the minimum 
of the effective potential. The size of the domain in this picture will depend on the time 
during which correlations had grown enough so that fluctuations start sampling beyond the 
spinodal point. 

Since this physical picture may not be studied within perturbation theory, we now intro- 
duce a non-perturbative method based on a self-consistent Hartree approximation, which is 
implemented as follows: in the initial Lagrangian write 

= ^ 2 (r,t))$ 2 (r-,t) + (~*V,*) - ^($ 2 (f,t))$ 2 (f,t)) (54) 
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the first term is absorbed in a shift of the mass term 



m 2 (t)->m 2 (t) + ^<$ 2 (t)> 

(where we used spatial translational invariance). The second term in (|53|) is taken as an in- 
teraction with the term ($ 2 (i))$ 2 (r, t) as a "mass counterterm" . The Hartree approximation 
consists of requiring that the one loop correction to the two point Green's functions must be 
cancelled by the "mass counterterm" . This leads to the self consistent set of equations 

/cl A* r d 1? 1 

j^y 3 (-iG£(M)) = J {2n)3 2 ^ <(A;) ^ + (tR-(t) coth[A^(fc)/2] (55) 



^ + P + m\t) + \^{t)) 



Ut = (56) 



Before proceeding any further, we must address the fact that the composite operator 
($ 2 (r, t)) needs one subtraction and multiplicative renormalization. As usual the subtraction 
is absorbed in a renormalization of the bare mass, and the multiplicative renormalization 
into a renormalization of the coupling constant. 

At this stage our justification for using this approximation is based on the fact that it 
provides a non-perturbative framework to sum an infinite series of Feynman diagrams of the 
cactus typef|, [36fl . 

It is clear that for t < there is a self-consistent solution to the Hartree equations with 
equation (|55|) and 



<$ 2 (t)> = <$ 2 ((T)> 

Ut = exp[ T ^<(*0] (57) 
A 

2 

(5* 



u\{k) = k 2 + m 2 + +($ 2 (0-)> = £ 2 + m 2 H 



where the composite operator has been absorbed in a renormalization of the initial mass, 
which is now parametrized as mf R = /^[(T 2 /T 2 ) — 1]. For t > we subtract the com- 
posite operator at t — absorbing the subtraction into a renormalization of m 2 which 
we now parametrize as rn 2 j R = — (T 2 /T 2 )]. We should point out that this choice of 
parametrization only represents a choice of the bare parameters, which can always be chosen 
to satisfy this condition. The logarithmic multiplicative divergence of the composite oper- 
ator will be absorbed in a coupling constant renormalization consistent with the Hartree 
approximation p6|, B7[, however, for the purpose of understanding the dynamics of growth 



of instabilities associated with the long-wavelength fluctuations, we will not need to specify 
this procedure. After this subtraction procedure, the Hartree equations read 

r d 3 k 1 

[<$ 2 (t)> - <$ 2 (0)>] = J N 3 2 , <(fc) ^(^W - 1] coth[/Wfc)/2] (59) 
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^ + fc 2 + m^) + ^(($ 2 (t))-($ 2 (0)) 



u£(t) = o 



m 2 R (t) 



T 2 

- 1 

J^2 



2 



1 } - 

J^2 



e(t) 



(60) 



(61) 



with T{ > T c and Tf <C T c . With the self-consistent solution and boundary condition for 
t < 



[<$ 2 (t < 0)) - <$ 2 (0)>] 
Z#(f <0) 





exp[=)=ic<j < (/c)t] 



/c 2 + m 2 



iR 



(62) 
(63) 

(64) 



This set of Hartree equations is extremely complicated to be solved exactly However 
it has the correct physics in it. Consider the equations for t > 0, at very early times, 
when (the renormalized) ($ 2 (t)) — ($ 2 (0)) ~ the mode functions are the same as in the 
zeroth order approximation, and the unstable modes grow exponentially. By computing 
the expression fl55|) self-consistently with these zero-order unstable modes, we see that the 
fluctuation operator begins to grow exponentially. 

As (($ 2 (t)) — ($ 2 (0))) grows larger, its contribution to the Hartree equation tends to 
balance the negative mass term, thus weakening the unstabilities, so that only longer wave- 
lengths can become unstable. Even for very weak coupling constants, the exponentially 
growing modes make the Hartree term in the equation of motion for the mode functions 
become large and compensate for the negative mass term. Thus when 



A 



R 



1m\ R 



($ 2 (t)> - <$ 2 (0)> 



the instabilities shut-off, this equality determines the "spinodal time" t s . The modes will 
still continue to grow further after this point because the time derivatives are fairly (expo- 
nentially) large, but eventually the growth will slow-down when fluctuations sample deep 
inside the stable region. 

After the subtraction, and multiplicative renormalization (absorbed in a coupling con- 
stant renormalization), the composite operator is finite. The stable mode functions will make 
a perturbative contribution to the fluctuation which will be always bounded in time. The 
most important contribution will be that of the unstable modes. These will grow exponen- 
tially at early times and their effect will dominate the dynamics of growth and formation 
of correlated domains. The full set of Hartree equations is extremely difficult to solve, even 
numerically, so we will restrict ourselves to account only for the unstable modes. From the 
above discussion it should be clear that these are the only relevant modes for the dynam- 
ics of formation and growth of domains, whereas the stable modes, will always contribute 
perturbatively for weak coupling after renormalization. 

Introducing the dimensionless ratios fl4*3| ) in terms of vtl^r ; m^R, (all momenta are now 
expressed in units of m/ ; ^), dividing (|60|) by rn 2 j R) using the high temperature approximation 
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coth[/3jo; < (A;)/2] ~ 2T i /u < {k) for the unstable modes, and expressing the critical temperature 
as T 2 = 24{i R /\ R , the set of Hartree equations ([59], |60|) become the following integro- 
differential equation for the mode functions for t > 



1+g dp 



with 



P 



2 ^jj[K{t)U p {t)-l] 



Kit) = o 



(65) 



U±{t<0) = exp[^<(g)t] 



R 



m 



2 \rp2 m2 

i,R _ L J j — l c 



m 



f,R 



[T 2 - Tf] 



4tt 2 



[T 2 _ T 2 



(66) 
(67) 

(68) 
(69) 



The effective coupling (^9|) reflects the enhancement of quantum fluctuations by high tem- 
perature effects; for Tf/T c ~ 0, and for couplings as weak as X R ~ 10~ 12 , g ~ 10~ 7 (Tj/T c ). 



The equations (J65|) may now be integrated numerically for the mode functions; once we 
find these, we can then compute the contribution of the unstable modes to the subtracted 
correlation function equivalent to ([47|) 



V^ HF \x,r) 



A 



R 



6m },R 



($(r,t)$(0,t))-($(f,0)$(0,0)) 



91 dp \7Ti? R , 



sin(px) 
(px) 



(70) 
(71) 



In figure (3) we show 



A 



R 



2m 2 fR 



(($ 2 (r)> - ($ 2 (0)» = 3(P"* (0, r) - X>"* (0, 0)) 



(solid line) and also for comparison, its zeroth-order counterpart 3(P^(0, r) — P^(0, 0)) 
(dashed line) for Ar = 10~ 12 , T,/T c = 2. (this value of the initial temperature does not have 
any particular physical significance and was chosen as a representative). We clearly see what 
we expected; whereas the zeroth order correlation grows indefinitely, the Hartree correlation 
function is bounded in time and oscillatory. At r « 10.52 , 3(D( hf \0,t) — T>( hf \0,t)) = 
1, fluctuations are sampling field configurations near the classical spinodal, fluctuations 
continue to grow, however, because the derivatives are still fairly large. However, after this 
time, the modes begin to probe the stable region in which there is no exponential growth. 
At this point A ff ((<fr 2 (r)) — $ 2 (0))), becomes small again because of the small coupling 



g « 10 , and the correction term becomes small. When it becomes smaller than one, the 
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instabilities set in again, modes begin to grow and the process repeats. This gives rise to an 
oscillatory behavior around A f ((<l )2 (r)) — $ 2 (0)) = 1 as shown in figure (3). In figure (4), 

f,R 

we show the structure factors as a function of x for r = 10, both for zero-order (tree level) 
(dashed lines) and Hartree V^ HF ^ (solid lines) for the same value of the parameters as 
for figure (3). These correlation functions clearly show the growth in amplitude and that 
the size of the region in which the fields are correlated increases with time. Clearly this 
region may be interpreted as a "domain", inside which the fields have strong correlations, 
and outside which the fields are uncorrelated. 

We see that up to the spinodal time t s pa 10.52 at which ^^((^C^)) ~ $2 ( ))) = 1> 
the zeroth order correlation 3P^(0, r) is very close to the Hartree result. In fact at r s , 
the difference is less than 15%. For these values of the coupling and initial temperature, 
the zeroth order correlation function leads to t s pa 10.15, and we may use the zeroth order 
correlations to provide an analytic estimate for t s , as well as the form of the correlation 
functions and the size of the domains. The fact that the zeroth-order correlation remains 
very close to the Hartree-corrected correlations up to times comparable to the spinodal is a 
consequence of the very small coupling. The stronger coupling makes the growth of domains 
much faster and the departure from tree-level correlations more dramatic [0. For strong 
couplings domains will form very rapidly and only grow to sizes of the order of the zero 
temperature correlation length. The phase transition will occur very rapidly, and our initial 
assumption of a rapid supercooling will be unjustified. This situation for strong couplings, 
of domains forming very rapidly to sizes of the order of the zero temperature correlation 
length is the picture presented by Mazenko and collaborators [O]. However, for very weak 
couplings (consistent with the bounds from density fluctuations), our results indicate that 
the phase transition will proceed very slowly, domains will grow for a long time and become 
fairly large, with a typical size several times the zero temperature correlation length. In a 
sense, this is a self consistent check of our initial assumptions on a rapid supercooling in the 
case of weak couplings. 

As we argued above, for very weak coupling we may use the tree level result to give an 
approximate bound to the correlation functions up to times close to the spinodal time using 
the result given by equation (|49"D, for Tf pa 0. Thus, we conclude that for large times, and 
very weakly coupled theories (Xr < 10~ 12 ) and for initial temperatures of the order of the 
critical temperature, the size of the domains £d(£) will grow typically in time as 



t D ® « (8^2)^(0)^/ ^ (72) 

with £(0) the zero temperature correlation length. The maximum size of a domain is ap- 
proximately determined by the time at which fluctuations begin probing the stable region, 
this is the spinodal time t s and the maximum size of the domains is approximately £o(t s ). 
An estimate for the spinodal time, is obtained from equation ( f49|) by the condition 
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3V(t s ) = 1. For weakly coupled theories and Tf ~ 0, we obtain 



t. 



y/2£(0) 



In 



'3A N 

47T 3 




(73) 



5 FRW Cosmologies: 

We now consider the case of scalar field in an homogeneous, isotropic and spatially flat FRW 
cosmology described by the metric 



ds 2 = dt 2 — a 2 {t)dx 2 
The action and Lagrangian density are given by 



S 

c 

K 



d 4 xC 



A 



a 3 (t) 

^[m 2 + ^]$ 2 (f,t) + ^$ 4 (f,t) 
\a cr / 



(74) 

(75) 
(76) 

(77) 
(78) 



with 1Z the Ricci scalar. We have introduced the coupling to the Ricci scalar as it will be 
induced by renormalization (see below). 

We can either follow the same method outlined above in terms of a path integral in the 
complex time plane, or alternatively one could evolve the density matrix in time by solving 
explicitly the Liouville equation, these two methods are equivalent. 

We have followed [|38| the second alternative in one-loop and Hartree approximation (thus 
providing another independent manner of tackling the problem). 

The initial condition is assumed to correspond to a local equilibrium thermal ensemble 
of the adiabatic modes at an early time t — > -co at a temperature T D = l//3 . In the 
case of a de Sitter background, this initial condition corresponds to a thermal distribution 
of the Bunch-Davies modes, and for T Q — > the density matrix describes a pure state, the 
Bunch-Davies vacuum. 

Any description of quantum statistical mechanics needs to specify an initial density ma- 
trix whose subsequent time evolution is determined by the Hamiltonian. The assumption 
of local thermal equilibrium at an early time is somewhat arbitrary. Whether or not it 
corresponds to a physically realistic situation can only be answered within the realm of a 
theory that incorporates gravity, particle physics and statistical mechanics as dynamical in- 
gredients. The choice of an initial condition for a density matrix will pervade any out of 
equilibrium dynamical calculation. If the system reaches thermal equilibrium after a long 
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time, the details on the initial conditions may ultimately be irrelevant, but this requires a 
deeper understanding. 

Our goal is to obtain the evolution equation for the order parameter, as well as the 
fluctuations. The order parameter is the expectation value (in the time dependent density 
matrix) of the volume average of the scalar field 



<j>{t) 



1 





d 3 x($(x, t)) 



1 

ft 



d 3 xTrp(t)$(x) 



cPafTrpiQU-^t, t )$(x, t a )U(t, t a ) 

(79) 

where Q is the comoving volume, and the scale factors cancel between the numerator (in 
the integral) and the denominator. U (t, t a ) is the time evolution operator with boundary 
condition U(t , t D ) = 1, t Q is the time at which the initial ensemble is specified. 



To one-loop order we find the effective evolution equations for the order parameter |38 

d*k put) +uut)\ 



+ 3^0 + n0) + A<4 



(2tt)3 2a*(t)W k (t 
where the mode functions U a ^ k (t) ; a = 1, 2 are real and satisfy 



coih\j3 o hW k (t o )/2] = (80) 



U ak — 

U\k{to) 



3 fa la., 

5 ; + ~IK*-t 
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3 a{t 



k 2 



a 2 (t) 



+ V"((j) cl (t))u ak = 







2a{t ) 



U 2k (t ) = W k (t ) 



a 2 (t ] 



+ V"(ct> cl (t )) 



(81) 
(82) 

(83) 



where 4> c i(t) is the solution to the classical equations of motion. 

We can see clearly that the time dependence of the mode functions is very important. If 
the scale factor varies rapidly with time as for example in a radiation dominated or de Sitter 
epoch, at no time will the notion of an effective potential be valid for the time evolution of 
the order parameter. 

The method lends itself to a self-consistent non-perturbative Hartree approximation. 
Invoking a Hartree factorization as in a previous section, we also derived the equations of 
motion. Introducing the fluctuation rj(x,t) = <&(x,t) — <p(t) and the Hartree "frequencies" 



V(%) = V''(0) + ^ 2 > 



we find 



+ 3-0 + V'(<P) + \<f>- 

(X — 

d 2 3 fa la 2 \ 



dt 



2 a 2 



h f d*k \u 2 k (t)+u 2 k (t)} 

(2tt)3 2a\t)W k (t ) 

+ v"m))+ 
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coth[/? aW fc (t )/2] = (85) 
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a 2 (t) 
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d*k [uut)+uut)\ 



(2tt) 3 
U lk {t ) = 1 



2a 3 (t)W k (t ) 
U 2k (t ) = 



coth [p KW h {t )/2] 



ak 
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U\ k (t ) 



3 a(t ) 



2 a(t ) 



Z4(t ) = W fc (t ) 



+ V (2) (0(t o )) 



The one-loop and Hartree equations present new features not present in our previous 
analysis in Minkowski space: first the "friction" term arising from the expansion, and sec- 
ondly the fact that physical wave-vectors are red-shifted. The red-shift will tend to enhance 
the instabilities as more wavelengths become unstable, but the presence of the horizon will 
ultimately constrain the final size of the correlated regions. 

We are currently studying the numerical evolutions of these equations and expect to 
report soon on details of the dynamics of the phase transition in these cosmologies [R3 . 



6 Renormalization Aspects 



To understand the renormalization procedure, it proves convenient to introduce the complex 
combination of mode functions 



u£(t) 



U lk =F iU 2k 
1 

3a(t Q 



2a(t ) 



(89) 
(90) 

(91) 



in the case of the one- loop approximation W k (4>) must be replaced by W k (<p). We will now 
analyze the renormalization aspects for the Hartree approximation, the one-loop case may 
be obtained easily from this more general case. We need to understand the divergences in 
the integral 

(27r) 3 2a 3 (t)W fc (t ) 1 ] 

The divergences in this integral will be determined from the large-k behavior of the mode 
functions that are solutions to the differential equation fl86|) with the boundary conditions 
( P0| , |9~Tp . The large-k behavior of these functions may be obtained in a WKB approximation 
by introducing the WKB function 



V k {t) 
V k {t ) 



exp 
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R{t')dt' 



satisfying the differential equation 



d 2 3 (a ld 2S 
M 2 ~~2 [a + 2^, 



^ r) +y"(m + ^(v 2 (x,t)} 



V k = 



(93) 
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with (eta 2 (x,t)) being the self-consistent integral in the Hartree equation (the one-loop 
approximation is obtained by setting this term to zero in the above equation). The mode 
functions of interest W m are obtained as linear combinations of the WKB function and its 
complex conjugate. The coefficients to be determined from the boundary condition at t Q . 
The function R(t) obeys a Riccati (WKB) equation 



R + R — 



a la 2 
a 2 a 2 







We propose a solution to this equation of the form 

R 



a(t) w k k 2 
and find the time dependent coefficient by comparing powers of k. We find 
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Finally we obtain 
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+ finite (101) 



where we have introduced a renormalization point K, and the finite part depends on time, 
temperature and K. For renormalization in the one-loop approximation, {rf) does not appear 
in the logarithmic divergent term in ( |101|) (as it does not appear in the differential equation 
for the mode functions up to one loop). There are several physically important features of 
the divergent structure obtained above. The first term (quadratically divergent) reflects the 
fact that the physical momentum cut-off is being red-shifted by the expansion. This term 
will not appear in dimensional regularization. 

Secondly, the logarithmic divergence contains a term that reflects the initial condition 
(the derivative of the expansion factor at the initial time t Q ). The initial condition breaks 
any remnant symmetry (for example in de Sitter space there is still invariance under the de 
Sitter group, but this is also broken by the initial condition at an arbitrary time t a ). Thus 
this term is not forbidden, and its appearance does not come as a surprise. 

The renormalization conditions in the Hartree approximation are obtained by requiring 
that the equation for the mode functions be finite|38[]. Thus we obtain 
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(103) 
(104) 



Notice that the conformal coupling £ = 1 /6 is a /ixed point under renormalization. 



7 High Temperature Limit: 

One of the payoffs of understanding the large-k behavior of the mode functions (as obtained 
in the previous section via the WKB method) is that it permits a straightforward evaluation 
of the high temperature limit. The finite temperature contribution is determined by the 
integral 

d*k [Uf k (t)+Ul(t)] 1 



//3 ° J (2tt) 3 a 3 (t)W k (t ) e^kito) _ i ( 105 ) 

For large temperature, only momenta k > T contribute. Thus the leading contribution 

is determined by the first term of the function R(t) of the previous section. After some 
straightforward algebra we find 
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Thus we see that the leading high temperature behavior reflects the physical red-shift, 
in the cosmological background and it results in an effective time dependent temperature 



Teff(t) = T 



a(t ) 
a(t) 



To this leading order, the expression obtained for the time dependent effective temper- 
ature corresponds to what would be obtained for an adiabatic (isentropic) expansion for 
blackbody-type radiation (massless relativistic particles) in the cosmological background. 

This behavior, however, only appears at leading order in the high temperature expansion. 
Higher order terms do not seem to have this form[^] as both the coupling to the curvature 
and the mass term distort the spectrum from the blackbody form. 



8 Conclusions: 

We provided a dynamical picture of the time evolution of a weakly coupled inflaton scalar 
field theory undergoing a typical second order phase transition in Minkowski space-time. 
The shortcomings of the usual approach based on the equilibrium effective potential and 
the necessity for a description out of equilibrium were pointed out. The non-equilibrium 
evolution equations for the order parameter were derived to one-loop and integrated for 
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"slow-rollover" initial conditions. We pointed out that the instabilities responsible for the 
onset of the phase transition and the process of domain formation, that is the growth of long- 
wavelength fluctuations, introduced dramatic corrections to the "slow-rollover" picture. The 
net effect of these instabilities is to slow even further the evolution of the order parameter 
(average of the scalar field) and for "slow rollover" initial conditions (the scalar field very 
near the false vacuum initially) the dynamics becomes non-perturbative. 

We introduced a self-consistent approach to study the dynamics out of equilibrium in this 
situation. It is found that for very weakly coupled theories domains growth to large sizes 
typically several times the zero temperature correlation length and that the growth obeys a 
scaling law at intermediate times. 

In strongly coupled theories the phase transition will occur much faster and domains will 
only grow to sizes of the order of the zero temperature correlation length. 

The one-loop equations of evolution for the scalar order parameter were obtained in 
spatially flat FRW cosmologies. These equations reveal that at no times is the approximation 
of an effective potential valid. We have also obtained the evolution equations in a non- 
perturbative self-consistent Hartree approximation. These equations present new features 
that illuminate the fact that the process of phase separation in cosmological settings is more 
subtle. 

The physical wave vectors undergo a red-shift and more wavelengths become unstable, 
thus enhancing the instability. On the other hand the presence of the horizon and the 
"friction" term in the equations of motion, will prevent domains from growing bigger than 
the horizon size. Eventually there is a competition of time scales that must be understood 
more deeply to obtain any meaningful conclusion about formation and growth of domains 



in FRW cosmologies. Work on these issues is in progress 38 



The renormalization aspects were studied and we point out that in a time dependent 
background, divergences appear that depend on the initial conditions. We also studied the 
leading behavior in the high temperature limit showing to this order that the temperature 
is red-shifted as in adiabatic expansion. 



This formalism could also be applied to chaotic initial conditions [p 
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